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We study three dimensional generalizations of the quantum spin Hall (QSH) effect. Unlike two 
dimensions, where the QSH effect is distinguished by a single Z2 topological invariant, in three 
dimensions there are 4 invariants distinguishing f 6 "topological insulator" phases. There are two 
general classes: weak (WT1) and strong (ST1) topological insulators. The WT1 states are equivalent 
to layered 2D QSH states, but are fragile because disorder continuously connects them to band 
insulators. The STI states are robust and have surface states that realize the 2+1 dimensional 
parity anomaly without fermion doubling, giving rise to a novel "topological metal" surface phase. 
We introduce a tight binding model which realizes both the WT1 and STI phases, and we discuss 
the relevance of this model to real three dimensional materials, including bismuth. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.20.-r, 85.75.-d 



In recent years, the advent of spintronics has motivated 
the study of the effects of spin orbit interactions (SOI) 
on the electronic structure of solids. SOI leads to the 
spin Hall effect 

which has been observed in GaAs|3, 
4]. We proposed [5| that in graphene the SOI leads to 
the quantum spin Hall (QSH) effect. In the QSH phase 
there is a bulk excitation gap along with gapless spin- 
filtered edge states. The QSH phase is distinguished from 
a band insulator by a Zi topological invariant 0, which 
is a generalization applicable to time reversal invariant 
systems of the TKKN invariant of the integer quantum 
Hall effect 0- Because of the weak SOI in carbon, the 
SOI induced energy gap in graphene is likely to be quite 
smallQ. However, Murakami has recently suggested that 
bismuth bilayers may provide an alternative venue for 
the QSH effect [9f. This breakthrough provides a new 
direction for the experimental observation of this phase. 

In this paper we consider the generalization of the QSH 
effect to three dimensions (3D). Our work builds on re- 
cent progress by Moore and Balents[l(|, who showed that 
time reversal invariant energy bands in 3D are charac- 
terized by four Z 2 invariants, leading to 16 classes of 
"topological insulators" . A similar result has also been 
obtained by Roy ll| . Here, we will explain the physical 
meaning of these invariants and characterize the phases 
they distinguish. One of the four invariants is of spe- 
cial significance and distinguishes what we will refer to 
as "weak" and "strong" topological insulators. With dis- 
order, the weak topological insulator (WTI) is equivalent 
to a band insulator, while the strong topological insula- 
tor (STI) remains robust. We show WTIs and STIs have 
surface states with an even and odd number of Dirac 
points respectively. The latter case leads to a new "topo- 
logical metal" surface phase, which we characterize. We 
introduce a tight binding model on a distorted diamond 
lattice, which realizes both the WTI and STI phases, al- 
lowing the surface states to be studied explicitly. 

In Ref. 12 we established the connection between the 



way which makes the generalization to 3D transparent. 
The 2D invariant can be understood using a Laughlin 
type construction [r| on a cylinder threaded by magnetic 
flux $ = or 7r (in units ofh/e). The invariant character- 
izes the change in the time reversal polarization (TRP), 
which signals the presence of a Kramers degeneracy at 
the ends, when $ is changed from to tt. If the cylin- 
der consists of a single unit cell in the circumferential (x) 
direction, then the magnetic flux threading the cylinder 
plays the role of the crystal momentum k x in band theory. 
The spectrum of the discrete end states of the cylinder 
as a function of flux then reflects the edge state spectrum 
as a function of momentum. The change in the TRP as 
a function of flux determines the way the Kramers de- 
generate end states at the edge time reversal invariant 
momenta (TRIM) k x = Ai = and k x = A2 = tt are 
connected to each other. In the QSH phase the Kramers 
pairs "switch partners" (Fig. la), reflecting the change 
in the TRP, while in the conventional insulator (Fig. lb) 
they do not. It follows that in the QSH phase edge states 
traverse the bulk energy gap, and cross the Fermi energy 
an odd number of times between A a and A;,. In the insu- 
lating phase, the edge states cross the Fermi energy an 
even number of times if at all. They are not topologically 
protected, since changes in the Hamiltonian at the edge 
can push the entire edge band out of the bulk gap. 

The change in the TRP between Ai and A2 is related 
to the bulk band structure, defined for a 2D system with 
periodic boundary conditions in both directions. The 2D 
Brillouin zone has four TRIM, Ti, which are related to 
— Ti by a reciprocal lattice vector. For an edge perpen- 
dicular to G the ID edge TRIM A a =i 2 are projections of 
pairs T i=a i, r i=a2 , which satisfy T al - T a2 = G/2, onto 
the line perpendicular to G. 

The TRP associated with A a can be expressed as TT a = 
Si = aiS i= a2, where[12| 



6i = V / Det[ w (r 4 )]/Pf[«;(r l )] = ±1. 



(1) 



Z 2 invariant for the bulk QSH phase and the spin filtered 
edge states. We begin by reviewing that argument in a 



Here the unitary matrix uiy(k) = k)|0|uj(k)). At 
k = Ti, Wij = —Wji, so the Pfaffian Pf [tt;] is defined. TT a is 
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FIG. 1: Schematic surface (or edge) state spectra as a function 
of momentum along a line connecting A a to A& for (a) ir a iTb = 
— 1 and (b) TT a TTb = +1. The shaded region shows the bulk 
states. In (a) the TRP changes between A a and As, while in 
(b) it does not. 



free of the ambiguity of the square root in (TT]), provided 
the square root is chosen continuously as a function of 
k. However, ir a is not gauge invariant. A k dependent 
gauge transformation can change the sign of any pair of 
Si's. This reflects the physical fact that the end Kramers 
degeneracy depends on how the crystal is terminated. It 
is similar to the ambiguity of the charge polarization [l2T| . 
The product, ttitt2 — 31626364, is gauge invariant, and 
characterizes the change in TRP due to changing the flux 
from Ai = to A2 = tt. This defines the single Z2 in- 
variant in 2D, and using the above argument, determines 
the connectivity of the edge state spectrum. 

In 3D there are 8 distinct TRIM, which are ex- 
pressed in terms of primitive reciprocal lattice vectors 
as r i= („ in2 „ 3 ) = (mbi + n 2 h 2 + n 3 h 3 )/2, with rij =0,1. 
They can be visualized as the vertices of a cube as in 
Fig. 2. A gauge transformation can change the signs of 
6i associated with any four Ti that lie in the same plane. 
Modulo these gauge transformations, there are 16 invari- 
ant configurations of <5j. These can be distinguished by 4 
Z2 indices vq', (y\i>2V 3 ), which we define as 

= n 5 «™ (2) 

nj=0,l 

(-1)— ,3 = FT 8nm _ (3) 

rij-^i— 0,l;ni — 1 

i>o is independent of the choice of b^. {v\L>2V 3 ) are not, 
but they can be identified with G u = ^ j; ^bi, which 
belongs to the 8 element mod 2 reciprocal lattice, in which 
vectors that differ by 2G are identified. (V1V2V3) can be 
interpreted as Miller indices for G„. 

^0-4 are equivalent to the four invariants introduced by 
Moore and Balents using general homotopy arguments. 
The power of the present approach is that it allows us 
to characterize the surface states on an arbitrary crystal 
face. Generalizing the Laughlin argument to 3D, con- 
sider a system with open ends in one direction and pe- 
riodic boundary conditions in the other two directions. 
This can be visualized as a torus with a finite thickness 
(a "Corbino donut"), which has an inside and an outside 
surface. Viewed as a ID system, we then seek to clas- 
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FIG. 2: Diagrams depicting four different phases indexed by 
va\ {y\v-2Vz). (a) depicts Si at the TRIM Fi at the vertices 
of the cube, (b) characterizes the 001 surface in each phase. 
The surface TRIM A a are denoted by open (filled) circles for 
7T a = S a i6 a2 = 1). They are projections of r a i and r a 2, 

which are connected by solid lines in (a). The thick lines 
and shaded regions in (b) indicate possible surface Fermi arcs 
which enclose specific A a . 



sify the changes in the Kramers degeneracy associated 
with the surfaces as a function of two fluxes threading 
the torus (or equivalently as a function of the two com- 
ponents of the surface crystal momentum). 

For a surface perpendicular to G, the surface Brillouin 
zone has four TRIM A a which are the projections of pairs 
Tai, Ta2, that differ by G/2, into the plane perpendicular 
to G. Due to Kramers' degeneracy, the surface spectrum 
has two dimensional Dirac points at A Q . The relative val- 
ues of 7r a = 6 a i6 a 2 determine how these Dirac points are 
connected to one another, as illustrated in Fig. 1. For 
any path connecting A a to A&, the surface band struc- 
ture will resemble Fig. la (lb) for 7r Q 7Tfc = — 1(+1), and 
the surface bands will intersect the Fermi energy an odd 
(even) number of times. It follows that the surface Fermi 
arc divides the surface Brillouin zone into two regions. 
The Dirac points at the TRIM A a with ir a = +1 are on 
one side, while those with 7r a = — 1 are on the other side. 

In Fig. [2] we depict 5i for four different topological 
classes, along with the predictions for the edge state spec- 
trum for a 001 face. The surface Fermi arc encloses either 
0(4), 1(3), or 2 Dirac points. When the number of Dirac 
points is not 0(4), there must be surface states which 
connect the bulk conduction and valence bands. 

There are two classes of phases depending on the par- 
ity of vq. For vq — each face has either 0(4) or 2 
enclosed Dirac points. For a face G = J^i m i^i there 
are 0(4) Dirac points for TOj = vi mod 2 (i = 1,2,3) 
and 2 Dirac points otherwise. These phases can be 
interpreted as layers of 2D QSH states stacked in the 
G„ direction. They resemble three dimensional quan- 
tum Hall phases [14j, which are indexed by a triad of 
Chern integers that define a reciprocal lattice vector G 
perpendicular to the layers and give the conductivity 
o~ij = (e 2 / 'h)eijkGk I '(27r) . In the present case, G^ is de- 
fined modulo 2G, so that layered QSH phases stacked 
along G„ and G^ + 2G are equivalent. 
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FIG. 3: Energy bands for (a) the model (g| with t = 1, Aso = 
.125. The symmetry points are T = (0,0,0), X = (1,0,0), 
W = (1,1/2, 0), K = (3/4, 3/4, 0) and L = (1/2, 1/2, 1/2) in 
units of 2ir/a. The dashed line shows the energy gap due to 
Sti — .4. (b) shows the phase diagram as a function of Sti 
and St2 (for bonds in the 111 and 111 directions) with phases 
indexed according to cubic Miller indices for G„ . The shaded 
region is the STI phase. 



The presence or absence of surface states in the i/q = 
phases is delicate. For the 0;(00T) phase in Fig. 2, the 
100 face has two Dirac points, while the 801 face has 
0(4). This sensitivity is a symptom of the fact that the 
topological distinction of these phases relies on the trans- 
lational symmetry of the lattice. Indeed, if the unit cell is 
doubled, the two Dirac points fold back on one another. 
A weak periodic potential then opens a gap. It is thus 
likely that disorder will eliminate the topological distinc- 
tion between these phases and simple insulators. Surface 
states will generically be localized. For this reason, we 
refer to the vo — phases as "weak" topological insu- 
lators. Nonetheless, the weak invariants have important 
implications for clean surfaces. 

The i/o = l phases are more robust, and we refer to 
them as "strong" topological insulators. In this case the 
surface Fermi arc encloses 1(3) Dirac points on all faces. 
If the Fermi energy is exactly at the Dirac point this pro- 
vides a time reversal invariant realization of the 2+1 di- 
mensional parity anomaly 15, 16|> 17, without fermion 
doubling. This can occur because the Dirac point part- 
ners reside on opposite surfaces. For a generic Fermi 
energy inside the bulk gap the surface Fermi arc will en- 
close a single Dirac point. This defines a two dimen- 
sional "topological metal" that is topologically protected 
because a quantized Berry's phase of n is acquired by an 
electron circling the Fermi arc. This Berry's phase im- 
plies that with disorder the surface is in the symplectic 
universality class TjJ 2(| , which is not localized by weak 
disorder. 

To develop an explicit model of these phases we con- 
sider a 4 band tight binding model of s states on a di- 
amond lattice with SOI, that generalizes the 2D honey- 
comb lattice model[l El [H|. 



The first term is a nearest neighbor hopping term con- 
necting the two fee sublattices of the diamond lattice. 
The second term connects second neighbors with a spin 
dependent amplitude, djj 2 are the two nearest neighbor 
bond vectors traversed between sites i and j. a is the 
cubic cell size. The energy bands are shown in Fig. 3a. 
Due to inversion symmetry, each band is doubly degener- 
ate. The conduction and valence bands meet at 3D Dirac 
points at the three inequivalent X points X r = 2irr/a, 
where r = X,y, z. This degeneracy is lifted by symmetry 
lowering modulations of the four nearest neighbor bonds 
t — ► t + St p , with p = 1, 4. 

Near X z the low energy effective mass model has the 
form of a 3+1 dimensional Dirac equation, 

H z cS = taa v q z + 4\ SO aa z (s x q x - s y q y ) + m z a x . (5) 

Here q = k X z and m z = J2 p <5i p sgn[d p • z], d p is the 
bond vector associated with the p'th nearest neighbor 
bond. The Pauli matrices a % are associated with the 
sublattice degree of freedom, while s l describe the spin. 

are ^ ne same with x, y and z permuted in qi and s*, 
but not a 1 . Transitions between different phases occur 
when the masses at any of the X r vanish. St p = is thus 
a multicritical point separating 8 different time reversal 
invariant phases. 

Determining Vi for these phases using (1-3) requires 
eigenvectors that are defined continuously throughout 
the Brillouin zone 12j . Since the Chern integers van- 
ish, this is always possible. Determining the appropriate 
phases numerically, however, is nontrivial. An alterna- 
tive (though tedious) numerical approach would be to 
characterize the Pfafhan function introduced in Ref. @. 
Generalizing the results of Ref. [13 it can be shown that 
the product of 4 Si in any plane is related to the ze- 
ros of the Pfaffian in that plane, which can be identified 
without choosing phases. For the present problem, how- 
ever, we are fortunate because the eigenvectors can be 
determined analytically, allowing for the continuation of 



, (k) between the different r,-. We find 



Si = sgn[E(* 



5t p ) cos Ti ■ (d p - di)]. 



(0) 



(ij) 



^E4s-(4 x 4) c ,- (4) 

({ij)) 



For small St p , 6 = 1 at k = and at 3 of the L points. 
5 = -1 at the 4th L point. At X r , 5 = sgn[m r ]. When 
one of the four bonds is weaker than the others (St p < 0, 
5t(p'/ p ) = for instance) the system is in a WTI phase, 
which may be interpreted as a QSH state layered in the 
dp direction. There are 4 such states, depending on p, of 
which two are shown in Fig. 3b. We labeled the phases 
with the conventional cubic Miller indices for G„. Ill 
and 111 are distinct elements of the fee mod 2 recipro- 
cal lattice. When one of the bonds is stronger than the 
others the system is in one of four STI phases. The band 
insulator 0; (000) is not perturbatively accessible from 
this critical point. However, in the tight binding model 
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FIG. 4: 2D band structures for a slab with a 111 face for the 
four phases in Fig. 3. The states crossing the bulk energy 
gap are localized at the surface. In the WTI (STI) phases 
there are an even (odd) number of Dirac points in the surface 
spectrum. The inset shows the surface Brillouin zone. 



it occurs when one bond is turned up so that 8t\ > 2t. 
A staggered sublattice potential also leads to a band in- 
sulator, but the strength must exceed a finite value (set 
by Aso) before that transition occurs. 

To study the surface states, we solve ((H) in a slab ge- 
ometry. Fig. 4 shows the 2 dimensional band structures 
of the four phases in Fig. 3 for a slab parallel to the 111 
surface along lines that visit each of the four TRIM. The 
plots with the same can be viewed as different faces of 
the same state. The bulk states above the bandgap are 
clearly seen. In addition, there are surface states which 
traverse the gap. In the WTI phases 0; (111) and 0; (111) 
there are and 2 2D Dirac points, on both the top and 
bottom surfaces, as expected from the general arguments 
given above. In the STI phases 1; (111) and 1; (111) there 
is 1(3) Dirac point on each surface. In each case, the non 
degenerate surface states near the Dirac points are spin 
filtered, such that (s(— k)) = — (s(k)). Spin density and 
charge current are thus coupled. 

Though the 4 band diamond lattice model is simple, 
it is probably not directly relevant to any specific mate- 
rial. However, it may give insight into the behavior of 
real crystals. Consider a sequence of crystal structures 
obtained from diamond by continuously displacing the 
fee sublattices in the (111) direction: 

diamond — > graphite(ABC) — > cubic. 

Starting with diamond, the 111 nearest neighbor bond 
is stretched, leading to the 0; (111) WTI phase. As the 
sublattice is displaced further both sublattices eventually 
reside in the same plane with a structure similar to ABC 
stacked graphite. Displacing further, the lattice eventu- 
ally becomes cubic. At this point, the gap closes, and 
the system is metallic. The s state model remains in the 
WTI phase up to the cubic point. 



Bismuth has the rhombohedral A7 structure, which 
can be viewed as a cubic lattice distorted "toward dia- 
mond" , along with a trigonal distortion of the fee Bra- 
vais lattice. Murakami showed that a bilayer of bismuth, 
whose structure is similar to a single plane of graphene, is 
in the QSH phase. This suggests that for weak coupling 
between bilayers bismuth is in the 0; (111) WTI phase. 
While this agrees with the simple model presented above, 
a realistic description of bismuth requires a theory which 
incorporates bismuth's five valence bands [2lj|. 

It will be interesting to search for materials in the STI 
phase, which occur on the "other side of diamond" in our 
sequence. We hope that the exotic surface properties pre- 
dicted for this phase will stimulate further experimental 
and theoretical efforts. 

It is a pleasure to thank Joel Moore and Leon Balcnts 
for helpful discussions. This work was supported by NSF 
grants DMR-0079909 and DMR-0605066 and DOE grant 
DE-FG02-ER-0145118. 

Note added: In subsequent work we have predicted 
that a number of specific materials are STI's[22j. These 
include the semiconducting alloy Bii-^Sb^ as well as a- 
Sn and HgTe under uniaxial strain. 
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